Quantum confinement of Mott electrons in ultrathin LaNiOs/LaAlOa superlattices 
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We investigate the electronic reconstruction in (LaNi03)„/(LaA103)3 (n =3, 5 and 10) superlattices due to 
the quantum confinement (QC) by d.c. transport and resonant soft x-ray absorption spectroscopy. In proximity 
to the QC limit, a Mott-type transition from an itinerant electron behavior to a localized state is observed. 
The system exhibits tendency towards charge-order during the transition, ah initio cluster calculations are in 
good agreement with the absorption spectra, indicating that the apical ligand hole density is highly suppressed 
resulting in a strong modification of the electronic structure. At the dimensional crossover cellular dynamical- 
mean-field calculations support the emergence of a Mott insulator ground state in the heterostructured ultra-thin 
slabofLaNiOs. 



PACS numbers: 

Prompted by the discovery of high-Tc superconductivity in 
cuprate compounds there has been a surge of activity to dis- 
cover materials with even higher transition temperature iT]. 
Recent remarkable advances in synthesis of artificial layers 
of complex oxides along with the progress in computational 
methods have re-energized the search for novel superconduc- 
tors outside of the cuprate family ||2] [3] . 

Towards the challenge, a recent proposal has been put for- 
ward to use heterostructuring and orbital engineering to turn 
hole-doped alternating unit-cell thin layers of a correlated 
metal LaNiOs (LNO) and band-gap dielectric LaAlOs (LAO) 
|4|. The proposal utilizes the Ni^^^ low-spin state in the 
bulk with a single unpaired electron occupying the degenerate 
(dx2_y2 /d^z2_,.2) eg orbital whose nodes point to the planar 
and apical ligands of the octahedra, respectively. Both t-J 
and LDAh-DMFT calculations [4 5 1 suggest that the quantum 
confinement together with the electronic correlations should 
make it possible to localize or empty the c?322_r2 band leaving 
the conduction electron in the dj.2_y2 band. Epitaxial strain is 
also suggested as a mean of the orbital control to manipulate 
the d^^2_^2 orbital to appear above d^2_y2 and play the anal- 
ogous role of the axial orbital of the high-Tc cuprates. 

Despite the attractive simplicity of the structure, the experi- 
mental realization of LNO/LAO superlattice (SL) presents a 
large degree of ambiguity in selecting a specific electronic 
ground state caused by (i) the intrinsic propensity of ortho- 
nickelates to charge-ordering ||6l |71, (ii) the orbital polar- 
ization due to chemical confinement imposed by interfacial 
bonding l|8]|9l, (iii) polar discontinuity ifTOlfTTl . (iv) epitaxial 
constraint on the crystal symmetry lfT2l [131 and (v) quantum 
confinement ||l4l[l5,|. Even for the undoped LNO/LAO SL 
(i.e. the proposed parent cuprate-like compound), the syner- 
getic action of these phenomena will likely render the the- 
oretical problem intractable to a priori predict whether the 
cuprate-like physics can be experimentally realized. 

Here we report on emergence of a Mott-type metal- 
insulator transition (MIT) at the dimensional crossover in the 



experimentally realized (LNO)„/(LAO)3 SLs. By using trans- 
port measurement and soft x-ray absorption spectroscopy in 
combination with cluster calculations and cellular dynamical- 
mean-field theory calculations, we demonstrate the presence 
of a charge-localized ground state at the quantum confinement 
limit caused by strongly enhanced correlations and broken 
translational symmetry across the interface. In sharp contrast 
to conventional bandwidth-control, the observed MIT is a 
continuous one with a critical region due to an emerging 
correlated gap imposed by confined dimensionality. The 
system also exhibits tendency towards charge-ordering as a 
competing state during the transition. Microscopically, the 
crossover from the localized to itinerant behaviors is driven 
by a strong competition between the reduced covalency of 
the Ni-O-Al bond and the recovery of the ligand hole density 
within the LNO slab. 

High-quality epitaxial [(LNO)„/(LAO)3]Ar SLs (n = 3, 
5 and 10 with = 10, 4 and 2 respectively) were grown 
by laser molecular beam epitaxy with in situ monitoring by 
Reflection High Energy Electron Diffraction, as described 
in Ref. lfTTl [T6ll : +2.2% tensile strain is applied with Ti02- 
terminated (001) SrTiOa substrates ifTTl . Atomic force mi- 
croscopy imaging reveals that the substrate surface morphol- 
ogy is well maintained after the growth; Symmetric scans and 
reciprocal space mapping by synchrotron-based x-ray diffrac- 
tion indicate coherently grown single crystal heterostructures 
with the fully strained state[16|. 

To explicate the evolution of the ground state due to quan- 
tum confinement, d.c. transport measurements from 400 K to 
2 K were performed in the van der Pauw geometry with a 
commercial physical properties measurement system (PPMS, 
Quantum Design). The obtained temperature dependent con- 
ductivity, cr(r), of the SL series is shown in Fig.[T|a). A tran- 
sition from metallic to insulating behavior can be clearly seen 
with decreasing n. The inset of Fig. [TJa) also demonstrates 
orders of magnitude change in a across the series taken at 
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FIG. 1 : (a) Conductivity of (LNO)„/(LAO)3 SLs versus tempera- 
ture. Inset: conductivity versus n at 10 K. The daslred line is a guide 
for eye. (b) Resistivity as a function of T^'^ at n = 10. The red 
solid line is a guid for eye. (c) dlno- as a function of XlT^I^ for 
the 71 = 3 SL shows the low temperature three dimensional variable 
range hopping behavior. 



10 K. As seen, in analogy to the bulk LNO, the n = 10 SL 
exhibits a highly metallic behavior ascribed to the three di- 
mensional regime. On the other hand, unlike the bulk which 
is known to follow a Fermi liquid behavior at low temper- 
atures ||18 I, the low-temperature resistivity p at n = 10 shows 
a peculiar c>c T^l"^ dependence below 200 K (see Fig. [IJb)). 
This power law dependence has been attributed to scattering 
by bond-length fluctuations |19| and indicates the proximity 
to a Mott-type metal-to-insulator transition (MIT) from the 
itinerant side due to enhanced electronic correlations in the 
quantum confinement regime. 

The intermediate thickness n = 5 SL still shows a clear 
metallic behavior at high temperatures but reaches a conduc- 
tivity maximum at T;nax= 120 K. While the metallicity of 
LNO films under large tensile strain is known to cease when 
the thickness is < 3nm [l20l], it is remarkable that the high 
temperature metallic behavior is well maintained in such a 
thin slab (< 2nm). This result unambiguously highlights the 
exquisite role of the confinement structure with inteifaces in 
defining the transport properties of SLs as compared to the 
bulk LNO and single layer films. Notice, the conductivity a 
here is also well below the loffe-Regel Umit El (~ 2000 
S/cm), manifesting the modification of the correlation nature 
of the observed behaviors and the minor role of disorder. Be- 
low Tmax, cr(?i„ax) — cr(2 K) (^381.6 S/cm) is as large as 
cr(2 K) and almost twice of a{T^^^ - (t(400 K). Neither 
hopping conductivity nor quantum corrections is found to be 
able to account for this da/dT > O region. This transitional 
behavior explicitly indicates that all the involved interactions 
and relevant energy scales, like electron coiTelations and the 
bandwidth, are of a similar strength and strongly competing 
with each other. Thus, rt = 5 falls into a critical region of the 
confinement-controlled MIT. Upon further decreasing n to 3, 
da/dT > at all temperatures and cr(T) is found to follow a 
variable range hopping behavior below 100 K (see Fig.[TJc)), 



indicating that the system enters the insulating phase at n = 3. 
The evolution with n shown above points to an opening gap 
in the excitation spectrum. 

To gain microscopic insight into the effect of confinement 
on the nature of the MIT, we carried out detailed soft x-ray 
absorption measurements. Since resonant x-ray absorption 
is local and element specific probe, it is well suited for in- 
vestigating electronic structure of buried layers. All spectra 
were acquired in the bulk-sensitive fluorescence yield mode 
at the Ni L-edge at the 4ID-C beamline of the Advanced Pho- 
ton Source f2!2\. Precise information on the Ni charge state 
was obtained by aligning all the spectra to a NiO standard 
which was measured simultaneously with the SLs. The en- 
ergy resolution was set at 100 meV. The normalized spectra 
obtained at 300 K are shown in Fig. |2|a,b). For direct com- 
paiison, the i-edge spectra of bulk LNO and SmNiOa |23 1 are 
included. As seen, the energy position and the lineshape of the 
SLs are close to those of the bulk RENiOa (RE = rare eaith) 
Il24ll25l . indicative of a Ni^^^ charge state of the 3rf^ low-spin 
configuration. On the other hand, while the multiplet effects 
are smeared in the metallic bulk LNO, the is-edge spectra 
of the SL show a pronounced two-peak structure (shaded in 
green), closely resembling the absorption spectrum of the bulk 
SmNiOa (see Fig. 2. (a)), which in the bulk is a well studied 
case for the charge-ordered (2Ni^^^^Ni^^H-Ni^^) insulating 
state ||6ll7]|. As n decreases this multiplet effect becomes more 
pronounced, particularly in moving from n = 10 to 5. The en- 
ergy separation of the two split peaks is also clearly increasing 
as shown in Fig. |2]^c). The coiTesponding trend at the L2- 
edge emerges as a developing low-energy shoulder at about 
870.5 eV (see Fig. 2.(b)). The observed evolution is remi- 
niscent of that reported for the bulk RENiOa when crossing 
the metal-insulator boundary 1251 . unambiguously revealing 
the carrier localization in the presence of a developing cor- 
related gap in the quantum confined LNO slabs and a latent 
tendency to charge ordering during the three-to-two dimen- 
sional crossover. In accordance with the observations from 
d.c. transport, the emergence of the multiplet even in the rela- 
tively thick LNO slab {n = 10) lends solid suppoit to strongly 
confinement-enhanced coiTelations. Meanwhile, when sub- 
ject to this enhancement, the inclination to charge ordering 
corroborates its role as a primary competing ground state not 
only in the bulk |6| but also in the predicted heterostructured 
LN0 14J. 

To shed additional light on the observed spectroscopic re- 
sults, we performed ab initio cluster calculations on a NiOg 
cluster. The model Hamiltonian includes on-site energies, 
crystal field. Coulomb interaction including the full multiplet 
structure and spin-orbit coupling |26|. The majority of the pa- 
rameters are equivalent to those in Ref |26|. The monopole 
pait of the Coulomb interaction is set at 6 eV. The major con- 
figurations used in the calculation are dg 



*-d„ I and d 



where d^^^ and L^gi stand respectively for a 3d hole and a lig- 
and hole with Cg symmetry and spin <j =^, I. Note that triva- 
lent nickel compounds, such as LNO, are significantly more 
covalent than divalent systems as the charge-transfer energy 



3 



CD 



o 



o 

< 




(b) 




\ 


















Vv n=-fQ 


















- ^ \ NBu/k 


















theoryyl 



> 

o 



o 
3 



03 



852 854 



856 870 

Energy (eV) 



872 874 




25 < 



3 4 5 6 7 
n (unit cells) 



10 



FIG. 2: X-ray absorption spectra of (LNO)„/(LAO)3 SLs at the La- 
edge (a) and L2-edge (b) for n = 3, 5, and 10 together with the 
theoretical calculations done for a NiOe cluster for a charge trans- 
fer energy of 0.8 (red) and 1.6 (blue) eV. Spectra of bulk LNO and 
SmNiOa 1 23] are included for comparison. The shaded (green) areas 
indicate the most significant changes during the evolution, (c) The 
observed multiplet splitting energy (square) and the corresponding 
calculated charge transfer energy (circle) versus n. Solid lines are 
fits to power-law. 



A„ = E{d"^^L) — E{d"-) decreases due to the change in 
electron count with Ay = Ag — J7, leading to strong mixing 
between the (f'L and rf^ configurations. The degree of Ni-0 
covalency here is varied by changing the charge-transfer en- 
ergy A7. Figure[2|a) shows the result of calculations when A7 
was varied from 0.8 eV to 1.6 eV. The calculation strongly 
implies that when the local Ni-O covalency is reduced with 
increasing A7, the observed evolution of the La-edge mul- 
tiplet splitting with n is well reproduced. To quantify the 
n-dependence, the value of A7 for each SL is obtained by 
matching the corresponding size of the splitting. As can be 
seen in Fig. |2]^c), the n-dependence of A7 clearly tracks that 
of the energy separation of the split peaks at the Ls-edge. 
With this excellent agreement between the experimental 



spectra and the calculation, we can speculate on a possible 
physical origin of the n— dependence of the correlated gap en- 
ergy. The increasing magnitude of A7 with decreasing num- 
ber of LNO unit cells can be connected to the effect of the 
confinement which breaks the translation symmetry with in- 
terfaces and causes the formation of a significantly less cova- 
lent Ni-O-Al chemical bond. The absence of d-electron states 
on Al and its much lower electronegativity compared to Ni 
strongly suppresses the ligand hole density on the apical oxy- 
gen Oa |8 9|. Consequently, because of the proximity to the 
more positively charged Al ion, the Madelung potential on 
the interfacial Oa is further raised with respect to the planar 
Op; this leads to further suppression on the apical hole density 
and charge transfer process. As a result, the interfacial Ni-Oa 
bond length becomes shorter H, the ligand hole is transferred 
into the Ni02 plane and the (PL configuration will acquire 
more of the dx2^y2 character This effect is in close analogy 
to the physics of cuprates where the Madelung energy of the 
apical oxygen and the Cu-Oq bond length are believed to be 
important parameters to control the planar hole density and to 
stabilize the Zhang-Rice states upon hole-doping lIZTl . 

As discussed before the cluster calculations can be success- 
fully applied to provide the insight how the local electronic 
structure is modified in proximity to the interface. However, 
the observed global MIT resulting from the competition be- 
tween the short-range electron correlation and the reduced 
bandwidth subject to quantum confinement is better treated 
within the layer extension of cellular dynamical-mean-field 
theory (layer CDMFT) [2 8 J with the Lanczos exact diago- 
nalization impurity solver 130 1. Note that CDMFT was re- 
cently applied to explain the opening of a correlation-induced 
pseudogap in SrVOs thin films 1.29 J . The layer CDMFT we 
employed here goes beyond that previous study by includ- 
ing short-range correlations within each monolayer to capture 
the physics of pseudogap. To make the problem tractable and 
yet relevant for the SLs we consider the multi-layered single- 
band Hubbard model for Ni 3d and simulate the high-energy 
Al 3s states as vacuum. The short-range correlation is anti- 
ferromagnetic as recently suggested in a similar system with 
the planar d^^-y^ orbital configuration |5 |. To neutralize the 
CDMFT's overemphasis on the in-plane correlations, we in- 
troduce a small anisotropy between the in-plane ti^ and out- 
of-plane tout transfers as ij„ = t{\ — 5) and to^t — ^(1 + 25) 
so that tout competes with the in-plane correlations. The 
bandwidth of the bulk non-interacting system is unchanged, 

w = m. 

Since bulk LNO is a metal we consider moderate values of 
U/t (U = 6t). Fig. [3] (a) presents the n-dependent evolu- 
tion of the spectral function. As clearly seen, the correlation 
gap in the two dimensional limit is filled up and becomes a 
pseudogap with increasing n, lending theoretical support to 
the notion that the experimentally observed MIT arises from 
the dimensional evolution of a quantum confined correlated 
carriers transforming from an insulator to a metal through a 
'bad' metal phase. Notice, the sharp quasiparticle peak is ab- 
sent at a; = even in the thickest slab because out-of-plane 



4 




-6 -4 -2 . 

0)/t 



1 2 3 4 5 e't'm 
n 



FIG. 3: Layer CDMFT results for the multilayer Hubbard model 
with various thickness n. The transfer anisotropy 5 is set at 0.2. (a) 
Average spectral function for U/t — 6. (b) Sheet conductance aa 
versus n for U/t = 5.5, 6 and 6.5. 



correlations are treated on a mean-field level. To see the re- 
lation between the (pseudo) gap opening and the anomalous 
transport properties, sheet conductance aa is computed by 
the standard Kubo formula. We introduce an imaginary part 
77 = O.lt to the real frequency cj to represent the finite impu- 
rity scattering. As shown in Fig.|3](b), with increasing n, an 
increases continuously and becomes practically unchanged 
above a critical thickness Uc below which the pseudogap is 
pronounced. These results reproduce well the critical regime 
at Uc w 5 in the transport measurements and strongly suggest 
the confinement-controlled short-range electron correlations 
as the origin of the observed MIT. As discussed earlier, it is 
known that bulk rare-earth nickelates have an intrinsic propen- 
sity towards the charge disproportionation when the electronic 
bandwidth is reduced [7J. The result is a complicated three- 
dimensional ordering of both charge and spin. While such an 
order might be thought as a complex competing state, we be- 
lieve that the essential physics of metal-insulator crossover by 
the quantum confinement is captured by the CDMFT. 

In conclusion, we have investigated the quantum confine- 
ment effects on the electronic properties of (LNO)„/(LAO)3 
SLs. By lowering the dimensionality of the heterostructured 
LNO slab, enhanced electron correlations are observed to 
drive the emergent Mott-type MIT with the latent competing 
state of charge-ordering. The strong multiple! splitting ap- 
pearing in X-ray absorption spectra is the effect of the inter- 
face leading to the redistribution of ligand hole density and 
reduction of Ni-O-Al covalency as confirmed by ab initio 
cluster calculations. The result indicates intriguing similari- 
ties between the physics of heterostructured nickelate-based 
Cg-system and parent cuprates. The critical region of the MIT 
at n = 5 deduced from the transport measurements is ascribed 
to the pseudogap opening due to dimensional crossover and 
strongly enhanced short-range correlations confirmed by the 
layer-resolved CDMFT calculations. 
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